Irrational vs. rational charge and statistics in two-dimensional quantum systems 
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We show that quasiparticle excitations with irrational charge and irrational exchange statistics exist in tight- 
biding systems described, in the continuum approximation, by the Dirac equation in (2+l)-dimensional space 
and time. These excitations can be deconfined at zero temperature, but when they are, the charge re-rationalizes 
to the value 1/2 and the exchange statistics to that of "quartons" (half-semions). 



Introduction - It was showii by Jackiw and Rebbi fl]] and 
by Su, Schrieffer, and Heeger I2i that excitations with fermion 
number 1/2 (or charge 1/2) exist at domain walls in the 
dimerization pattern of electrons hopping along a chain, as 
is believed to occur in polyacetylene. For electrons hopping 
on the honeycomb lattice, the lattice relevant to graphene, a 
topological defect in a dimerization pattern that is realized by 
a vortex was shown to lead to a topological zero-mode and 
bind the fermion number 1/2 to the vortex fstj- Fermion- 
number fractionalization in both polyacetylene and graphene 
can be understood in terms of the spectral properties of one- 
dimensional (ID) and two-dimensional (2D) massive Dirac 
Hamiltonians, respectively, that describe the low energy limit 
of the electronic tight-binding Hamiltonians. These fraction- 
ally charged topological excitations are generically decon- 
fined in ID. Their deconfinement in 2D relies on a mechanism 
for the screening of the 2D Coulomb potential by thermal fstl 
or quantum fluctuations involving an axial gauge field 

Applying different potentials to odd and even sites of 
the linear chain results in a continuously varying fractional 
fermion number S |^ 01 ■ At the level of the Dirac Hamil- 
tonian, this perturbation is represented by a second (gap- 
opening) mass term, that adds in quadrature to the mass due 
to the hopping dimerization of the chain. A complex order 
parameter is constructed from these two masses as its real and 
imaginary pieces, and the fractional fermion number is related 
to the phase twist of this order parameter as it sweeps through 
a domain wall, a result that has a natural interpretation within 
a bosonization scheme jst]. Although exchange statistics is 
ill-defined in ID, this varying phase twist also implies a con- 
tinuously varying exclusion statistics 

Can the charge and the exchange statistics of fractionalized 
quasiparticles in 2D be continuously varied as well? Here 
we show that they can. The fermion numbers of quasiparti- 
cles bound to a vortex can thus be irrational. Remarkably, 
if an axial gauge field supporting a half vortex is added to 
(precisely) screen the interaction potential between quasipar- 
ticles, their fractional fermion number re- rationalizes to the 
value Q = 1/2 and their statistical angle (for when time- 
reversal symmetry is broken) to the value 9/tt = 1/4. These 
results are first derived at the level of an effective field theory 
in (2 + 1) -dimensional space and time. We then discuss the 
relevance of this analysis for planar tight-binding models. 



Definitions - The massive Dirac Lagrangian that we shall 
consider in this letter takes the form 

£ [r {^^^ + 75^5m) - ^, (la) 

with /i = 0, 1, 2 and a = 0, 1, 2, 3, ^ = 1-^7°, and 4 x 4 
matrices 
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We allow for the background fields A^^ and to vary in 
space-time x = (x^) = (t,r). The field couples to the 
Dirac fermions as an axial U{1) gauge field does 10]. The 
four fields A^, when constant in space and time, open an 
energy gap in the Dirac spectrum. These masses have the 
following physical meaning on the honeycomb lattice, e.g., 
graphene. The masses Ai 2 coiTespond to the two compo- 
nents of the complex Kekule bond-density-wave order param- 
eter A = A]^ + zA2 jst). The mass A3 = is a staggered 
chemical potential that favors charges to sit in one of the sub- 
lattices of the honeycomb lattice ||9[]. Finally, Aq = ?] is the 
only of the four masses that breaks time-reversal symmetry 
(TRS). It originates on the lattice from a next-neighbor hop- 
ping term with phases that was introduced in Ref. We 
will use a constant ?/ to discuss quasiparticle statistics. The 
TRS masses fi^, A^, and Aj add in quadrature. Hence, we 
define the fields m and n = {rii, rig) by 

m := jAf + Aj+fil n, := ^. (2) 
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Induced fermionic U{\)-charge current with TRS - The 
current = {^{x)^^'^{x)) can be computed perturba- 

tively in a derivative expansion in n. When = 0, it reads 

f = l:;^e^^'"'n■{^,n^^^n). (3) 

Equation (O is the desired extension to 2D of the bosonized 
cuiTent in a Luttinger liquid. For a static charge- 1 vortex that 
vanishes at the origin and has the large distance asymptotic 
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FIG. 1: (color online). The 2 sphere spanned by the masses 
n = (ni, 712,713) — m^^(Ai, A2, ^s). A vortex of charge 1 in 
the dimerizations Aj and A2 corresponds to a parallel on the sphere 
whose height is fixed by the asymptotic value of the staggered chem- 
ical potential (1^(00). The spherical area red (or yellow) enclosed by 
this parallel gives the fractional charge induced by the fermionic zero 
mode localized at the core of the vortex. 



form 



A{r) = A{oo)e"f' + Oir-^), 



(4) 



where r = r(cos (j>, sin (j>), we find the charge 
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Ms(oo)/^(00) 

Vl + M.2(oo)/A2(oo) 



(mod 1) (5) 
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that varies continuously with the asymptotic value of the di- 
mensionless ratio /ig(oo)/A(oo). The ambiguity mod 1 in the 
fermion number Q corresponds to whether the bound state in 
the gap is empty or occupied. For yUs(oo) = 0, (5 = ±1/2 is 
recovered as in Ref . [30 . 

^ 0, the induced fermionic cuiTent is = 
Sjf^ where j^^^ is Eq. (O in which the derivatives have 
been replaced by the covariant derivatives 13^ ~ + '2iA'^, 
Sj'' = - 1 ajjjj ^ e''''Pd„A^,. The current > 

is conserved. The additional charge contributed by Sj^ is 
proportional to the flux carried by Ag^. In fact, in the pres- 
ence of the static vortex ^ and of the static half-vortex 
A|(r) = -t'^^a^{r) where 05(0) = and 05(00) = 1/2, 
the fractional charge is Q = 1/2 for any /ig(cxD)/A(oo)! 

Observe that, when Ar^^ = 0, despite the fact that the 
Dirac Hamiltonian associated to Eq. ([T]l does not possess an 
SU{2) symmetry (the mass-matrices Mi 2 3 do not close an 
SU{2) algebra under multiplication), the induced fermionic 
current ^ is the same as that in the 0(3) nonlinear sigma 
model (NLSM) derived from Dirac Hamiltonians with an in- 

This is so for the fol- 



m 



ternal SU{2) symmetry yj. 
lowing reason. The nonunitary transformation ^ = x^'Iq 
and = X is harmless as it induces a constant Jacobian 
that cancels for any physical observable. Under this trans- 
formation, the matrices = MqAI^ close the SU{2) alge- 
bra [Sj, E^] = 2iej^j.E^ and commute with the transformed 
matrices = -^^o7^- The transformed Lagrangian ([T]i is 
then invariant under the SU (2) transformations generated by 
S = (El, E2, S3) when A'^ ^ 0. 

Exchange statistics - If two identical quasiparticles associ- 
ated with two vortices are exchanged, the many-particle wave- 



function changes by the phase cxp(z6') with 6 the statisti- 
cal angle. On general grounds, we expect that fractionally 
charged quasiparticles acquire a mutual fractional statistics. 
However, for quasiparticles to display fractional statistics one 
needs to break TRS, and thus we now include a uniform TRS- 
breaking Aq = The calculation of the fractional charge 
when 77 = has shown that fractionalization is of topolog- 
ical origin in that it only depends on the asymptotic values 
of the fields. The same fractional charge follows when we 



impose the nonlinear constraint n = 



= 1, 



with 777 in Eq. (|2]i held constant (see Fig.[T]i. Now, a uniform 
77 competes with a uniform 777 leading to a quantum critical 
point when m = rj. Correspondingly, the fractional charge 
can be computed for a nonvanishing ratio 77/777 and shown 
to be the same as before except for the multiplicative factor 
8(777 — 1 77 1) with O the Heaviside step function. More- 
over, for A^ = 0, the effective action governing the dynam- 
ics of the background fields n is, despite the lack of SU (2) 
symmetry in Eq. ([T]i, the same as that of the 0(3) NLSM of 
Refs. lfTll[l2ll: 



cff 



8tt 



^ J d^x {d^nf + inn^ sgii(77) S'nopf + • • • (6a) 



The Hopf term is nonlocal in terms of the fields n. One pos- 
sible representation is 11311 (see also Ref. 1 1 4) 



d^xtv [{Uid'^U'') {Uid-'U'') {UiOm'')] 

(6b) 

with U the unitary 4x4 matrix-field that rotates the space- 
time dependent vector n into a fixed direction, say (0, 0, 1): 
H ■ n — U'^J^^U. When the angle Kjj :— 8(|77| — 777) is 
nonvanishing, the angle vanishes and vice versa. Both the 
fractional charge and the Hopf term contribute to the statistical 
angle. 



^Hopf 



e/TT -=sgn(7/) Q^(^Kc + = sgn(77) 



(7) 



A heuristic derivation of Eq. ^ when 777 > jr/j and A^ ~ 
is the following. Consider first a uniform n, for which the 
square of the Dirac Hamiltonian H associated to Eq. ^ is 



2riJ2i=i '^i^i- This Hamiltonian 



resembles that of a particle with spin in a magnetic field, of 
strength proportional to 77, that points in the n direction (here- 
after, we take 77 > 0). If the uniform n is rotated adiabati- 
cally by 27r around the z direction, a Berry phase 7^.^^ (n) = 
27r sin^ ^ is accumulated for each filled single-particle elec- 
tronic state, where a is the angle that n makes with the north 
pole of the unit sphere. The many-body Berry phase is the 
sum over occupied single-particle Berry phases after subtrac- 
tion of the Berry phase of some reference many-body state 
(i.e., without vortex). Now, n is not spatially uniform for a 
vortex, and thus if all n{r) are rigidly rotated by 27r, their ac- 
cumulated phase depends on the texture. In a semiclassical 
approximation, the Berry phase at each r must be weighted 
by the fermion density: 7j.qj = / rf^r j°(r) 27rsin^ 
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FIG. 2: (color online), (a) The columnar dimerization pattern, as 
indicated by the coloring of the bonds, for the nearest-neighbor hop- 
ping on a square lattice in the background of the yr-flux phase that 
opens up the gap A(oo). A staggered chemical potential, as indi- 
cated by the distinction between sites of sublattice A and B opens 
up the gap ^g(co). (b) The staggered dimerization pattern for the 
nearest-neighbor hopping on a square lattice that generates an axial 
gauge field, (c) A charge-1 vortex in the columnar dimerization (a), 
(d) A charge- 1/2 vortex in the staggered dimerization (b) with core 
radius c = 0.3. (e) The fermion density profile of (c) for a square 
lattice with open boundary conditions and 144 x 144 sites, (f) The 
fermion number as a function of the scaling variable /is(oo)/A(co) 
in the presence of the single charge-1 vortex (c) or with the addition 
of the axial charge-1/2 vortex (d) with core radius c = 0.01. The 
staggered chemical potential jj,^ takes the values O.Olt (black), O.OSt 
(red), 0.06t (green), and O.lt (blue). The thick and thin lines are the 
prediction lO without the axial vortex and Q — 1 /2, respectively. 



With the help of Eq. (O when m > rj, this gives 7^.^^ = 
/ IP 27rsin^ f = Trsin"* f = ttQ^. Hence, the phase accu- 
mulated by spinning the vortex by 2tt is tt Q^, which is what 
an object with spin collects. The exchange statistics for 
two quasiparticles of charge Q each should then be ^/tt = Q^. 

Numerics - We compare the predictions that the charge 
varies continuously as a function of the scaling variable 
/ig (00) /A(oo) according to Eq. ^ in the absence of the axial 
vortex [05(0x0) = 0] and that the charge re-rationalizes to Q = 
1/2 in the presence of the axial half-vortex [05(00) — 1/2] to 
results from an exact diagonalization in a lattice model. The 



model we consider is the 7r-flux phase lIlSll for electrons hop- 
ping on a square lattice, to which we add a dimerization pat- 
tern of the hopping amplitudes that realizes, on the lattice, the 
mass and axial vortices of Refs. Qj. We then study the 
charge bound to these singularities as a function of the stag- 
gered chemical potential /.t^ for the cases 05(0x0) = 0,1/2 (see 
also ifiell). 

Consider a square lattice as in Fig. |2] on which spin- 
less fermions hop. The square lattice can be divided 
into two interpenetrating sublattices A (open circles in 
Fig. |2]i and B (filled circles in Fig. |2]i, such that all near- 
est neighbors of sites in A belong to B, and vice versa. 
Here we construct the tight-binding Hamiltonian as H = 



+ H.C., where r 



(7711,771,2) G labels the sites in A, located at positions 
{nil + ™2)^ + ("^1 ^ '^2)y ill square lattice. The 
are the four vectors (labeled counterclockwise, starting from 
the +x direction, hy j = 1, 2, 3, 4) connecting a site in A to 
one of its four nearest-neighbor sites in B. 

The kinetic energy in Eq. ([T]) follows from linearizing at 
half-filling the tight-binding dispersion when the only non- 
vanishing hopping amplitudes in units of t > are f!^ j = 
(_l)("h+"2)(*,,i+*,,3)t. These hoppings define the tt-Aux 
phase (they are gauge equivalent to the case of uniform hop- 
pings t but with a magnetic flux of tt in units of h/ e threading 
each elementary plaquette, which is indicated in Fig.|2]). 

A staggered chemical potential on sublattice A and 
— /ij, on sublattice B induces the bilinear ^M^'^ in Eq. ([T]l |@]. 
A dimerization, such as shown in Fig.|2|a), arises from 5t^ j cx 
(-l)™i+'n2(^_^. ^ _ 6j^^)tl,y Such dimerizations induce the 
mass bilinears 2^ in Eq. ([T]) [3]. A dimerization, such 
as shown in Fig.l^b), arises from 5t^ ,j oc ((5^- 3 — Sj i)t'^.j. 
Such dimerizations induce the axial bilinears 275^ i" 
Eq. O. 

A charge-7i vortex in the columnar dimerization pattern is 
defined by 

St,^ „ = C ., — (cos n0 cos (f>j — sin 776* sin 0^ ) , 



cos 6 



2t 



(8) 



- = 7r(777i + 7772) + j(7r/2), 



and shown in Fig.|2|c) when 7i = 1. A half vortex with core 
radius c in the staggered dimerization pattern is defined by 



5t 



r,3 



'''-,7 

2\r\ 



tanh — [(5^- 3 - S^j) cos0 + (5^ ^ - 8^ 2) sin^] 

(9) 

and shown in Fig.|2jd). 

The charge-density profile is obtained by adding the con- 
tributions from all exact single-particle eigenstates that are 
filled: negative energy states plus the bound state with en- 
ergy within the gap that appears as a consequence of the vor- 
tex. The fermionic density profile of the vortex from Fig.|2c) 
is shown in Fig. [2|e) for a lattice made of 144 x 144 sites 
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with a weight concentrated either at the core of the vortex or 
on the boundary (open boundary conditions are used). The 
fermion number in the presence of a vortex is approximated 
by integrating the local fermion density in a disk surround- 
ing the vortex that extends beyond the localization length of 
the induced bound state but remains insensitive to the fermion 
density that has accumulated at the boundary. Subtraction 
of the background fermion charge in the absence of the vor- 
tex is always implicit. The continuous dependence of the 
fermion number ^ on the scaling variable /ig(oo)/A(oo) is 
shown in Fig. |2f) for the charge 1 vortex of Fig. |2c). Also 
shown in Fig. |2f) is the effect of the axial charge- 1/2 vor- 
tex from Fig. |2d) superimposed to the charge 1 vortex of 
Fig.|2c); as anticipated we find that the fermion number stays 
close to 1/2. We conclude that the agreement between the 
field theory and the numerics is good and improves as both 
//g(cx))/A(oo) and A(cx))/i -> 0. 

Energetics - At the level of the field theory, the energy cost 
to separate vortices of opposite charge grows like the loga- 
rithm of their separation. Temperature fluctuations screen this 
interaction through the Kosterlitz-Thouless mechanism. Al- 
ternatively, it was shown in Ref. 10] that the bare interaction 
between vortices is screened at zero temperature by coupling 
them to an axial gauge field A^^ carrying a half-vortex. The 
issue of energetics on the lattice is more subtle than in the con- 
tinuum approximation. The axial U{1) symmetry of the field 
theory ([T]i is reduced to the finite group Z4 of rotations by 
tt/2 about a site of the square lattice when the dimerization 
pattern responsible for the gap at the Fermi energy is com- 
mensurate with the lattice. If so, the energy cost to create two 
vortices of opposite charge grows linearly with their separa- 
tion. However, slight deformations of the hopping amplitudes 
away from the Z4-symmetric ones move the position of the 
Dirac points in the Brillouin zone, rendering the wavevector 
for the gap-opening dimerization pattern (which connects the 
two Dirac points) incommensurate with the lattice while pre- 
serving the bound state in the gap [Isl]. This incommensura- 
tion can restore an energy cost for the creation of a vortex- 
antivortex pair with a logarithmic dependence on their separa- 
tion, as is believed to happen in quantum dimer models, in the 
mechanism termed Cantor deconfinment in Ref. iflTil . 

Summary - We have calculated the irrational fermion num- 
ber attached to vortices in a complex-valued Higgs field cou- 
pled to massive Dirac fermions in (2 + 1) -dimensional space 



time. We also calculated the exchange statistics of these vor- 
tices and showed that it is irrational as soon as the fractional 
fermion number is. Deconfinement of the vortices at zero 
temperature requires the (axial gauge) coupling to a half vor- 
tex IJB . Remarkably, this coupling re-rationalizes the fermion 
number to the value 1/2. We have compared our predictions 
for the induced fermion number with an exact diagonaliza- 
tion study of a planar tight-binding model with a static defec- 
tive dimerization of the hopping amplitude and found a good 
agreement. 
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